Fluctuation Conductivity in Insulator-Superconductor Transitions with Dissipation 
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We analyze the fluctuation conductivity near the critical point in a 2D Josephson junction array 
shunted by an Ohmic resistor. We flnd that at the Gaussian level, the conductivity acquires a log- 
arithmic dependence on T/m (m, the inverse correlation length) when the dissipation is sufficiently 
weak. In the renormalized classical regime, this logarithmic dependence gives rise to a leveling of the 
resistivity at low to intermediate temperatures when fluctuations are included. We show, however, that 
this trend does not persist to T = at which point the resistivity vanishes. 



In two dimensional superconductors, the resistivity 
vanishes not at the mean-field temperature, TbcS: 
which the pair amplitude is established but at a lower 
temperature, Tc, where global phase coherence obtains 
[01 . The prominence of phase fluctuations in 2D is pri- 
marily responsible for the discrepancy between the on- 
set of pairing and the subsequent thermodynamic phase 
transition to a state with zero resistance. In 2D the zero- 
resistance state posseses algebraic long-range order . 
Above Tc but below TbcSj the system is in a paraco- 
herent state characterised by incoherent motion of the 
Cooper pairs. In this regime, phase fluctuations give rise 
to the conductivity which at the mean-field level is of the 
Aslamazov-Larkin form jsj . 

We are concerned in this work with the role of dissipa- 
tion on the conductivity in the regime where phase fluc- 
tuations dominate. While thermally-excited vortices are 
the most common dissipative mechanism in 2D supercon- 
ductors in the presence of a magnetic field, we specialize 
to the zero field case and consider here the case of shunt- 
ing an array of Josephson junctions with an Ohmic re- 
sistor Resistively-shunted Josephson junction arrays 
P-|l4| in 2D represent an idealization of an important 
class of insulator-superconductor systems, namely granu- 
lar superconductors [p^-pTf in which the role of Coulomb 
interactions and dissipation is paramount. At T = 0, 
such systems undergo an insulator-superconductor tran- 
sition (1ST) when the Josephson coupling energy exceeds 
the Coulomb charging energy. Dissipation suppresses 
quantum fluctuations and hence favours the supercon- 
ducting state. However, in the paracoherent regime, dis- 
sipation can also slow the motion of Cooper pairs. Naive 
considerations suggest that the fluctuation conductivity 
should decrease by an amount proportional to the inverse 
of the strength of the dissipation. In fact. Mason and Ka- 
pitulnik have recently suggested that dissipation in 
the presence of disorder can give rise to a leveling-off of 
the resistivity at T = 0-that is, a "metallic state" . 

We calculate here the role dissipation and fluctuation 
effects play in the fluctuation conductivity in 2D Joseph- 
son junction arrays. Within a Gaussian theory, we show 



that the standard [g form for the fluctuation conductiv- 
ity does not hold when the strength of the dissipation is 
less than the gap proportional to the inverse correlation 
length. In this regime, wc find that quite generally the 
conductivity scales logarithmically with the strength of 
the gap. While calculations of the fluctuation conduc- 
tivity have been performed previously none have re- 
vealed the regime with the logarithmic behavior at finite 
temperature. We then analyze the effect of the quar- 
tic term on the fluctuation conductivity in the N = oo 
limit. In a strictly D = 2 system, the N ^ oo limit 
results in ordering only at T = jl^. We show that in 
the renormalized classical regime, the logarithmic tem- 
perature dependence is transformed into a constant re- 
sistivity at low to intermediate temperatures. However, 
at sufficiently low temperatures, the resistivity is shown 
to vanish exponentially. 

To formulate the conductivity, we replace the micro- 
scopic Hamiltonian for an array of Josephson junctions 
with an effective Ginsburg-Landau theory. The coarse- 
graining approximation of Doniach offers a straight- 
forward way of obtaining the appropriate Landau the- 
ory. To this end, we introduce the complex order pa- 
rameter ?/'(r, r) whose expectation value is proportional 
to (exp(i0)), where (f> is the phase of a particular junc- 
tion. An effective Landau theory is valid if ip is suffi- 
ciently small as is the case in the vicinity of the onset of 
global phase coherence. The minimal Ginsburg-Landau 
theory [f7|- p^ , pO[ required to model quantum fiuctuations 
and dissipation near the zero-resistance quantum critical 
point is the Gaussian free energy functional. 
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where A(r, r) is the vector potential, e* = 2e, is 
proportional to the inverse correlation length, and k and 
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r] measure the strength of the quantum fluctuations and 
dissipation, respectively. In Fourier space, the dissipation 
term, Ldis = vJ2k,uj„ l^n||V'(k, w„)p, corresponds to the 
phenomenological Ohmic model introduced by Caldeira 
and Leggett |^. We have retained the vector potential, 
A(r,T), to facilitate the calculation of the conductivity. 
Along the imaginary frequency axis, the conductivity 
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is related to the fourier transform of the second variation 
of the partition function Z with respect to the vector 
potential. The partition function is defined in the usual 
way as the functional integral of exp(—_F [■)/)]). 

In the zero wavevector limit, the conductivity to 1-loop 
order [E0| can be rewritten as 



(G(k,c^„)-G(k,u;, 



2klGik,u;.ra) 
. + w„)) , 



(2) 



where G(k, lj„) is the standard Matsubara Green func- 
tion. To obtain the conductivity for real frequencies, 
we must analytically continue to real frequencies. As 
van Otterlo, et. al. have pointed out, inclusion 
of dissipation changes the analytical properties of the 
Matsubara sums. We define the retarded and advanced 
Green functions, G^{z) = {k'^ + m? — k^z^ — ij]z)~^ and 
G^{z) — (fc^ + — K^z^ + irjz)^^ . Upon choosing the 
appropriate contour, we find that 

^{^) ^y^l k'dk J ^ coth —dz [(G^(z) - G^(z)) 

[G^iz) + G^{z)-G^iz + uj)-G^iz~u)]] . (3) 

We are particularly interested in the limit of zero fre- 
quency. In this limit, the product of Green functions 
simplifies significantly leading to 



a{Lo = 0) = 
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k^dk 



dx 



sinh X 



[(e|-4T2K2a;2)2+ 42^2^22-2] 



(4) 



as our working expression for the zero-frequency conduc- 
tivity. In this expression, e\ = k"^ -\- m? . 

From the analytical structure of the integrand in 
Eq. (Q), it is clear that distinctly different regimes arise 
when 1) m <C nT and 2) ra':^ kT. In each of these cases, 
we explore the limiting form for the conductivity as the 
magnitude of the dissipation varies. In our treatment, we 
explicitly assume that ti/k <^ 1. In the limit of 77 = 0, the 
pole in Eq. makes the conductivity infinite at least to 
1-loop order as obtained previously by van Otterlo, et. 
al. (20| and Damle and Sachdev pl[ . Two loop and higher 



order corrections regularize the static conductivity even 
when 7] = 0, however. Nonetheless, as the dissipation 
is turned on, the divergent 1-loop conductivity becomes 
finite. In the present treatment, we address only the 
form of the static conductivity at lowest order. We focus 
first on the region characterized by to ^ kT. The two 
subregimes of interest are a) rj/n <^ m and b) m <C tj/k. 
Only the latter regime has been considered previously [|j . 
It is in the former regime that the new logarithmic be- 
haviour emerges. To obtain this result, we note that the 
main contribution from the integral in Eq. arises from 
small X and the minima of the denominator of the inte- 
grand which occur at x| = l/4K'^T'^[el - /2T'^]. The 
latter two roots are real only if to > ti/k. If we work 
in the limit in which to 3> rj/ k,^ we expand the integral 
about x± and perform the x integration 
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ed^ 
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by introducing the new variables, ^ — k/2K.T and m = 
to/2kT. In the limit that to ^ 1, we can approximate 
sinh^ ^ + m? as ^2 ^ 777,2 ^ jg approximation is valid 
if <C 1 or equivalently an ultraviolet cutoff is intro- 
duced on the momentum of order 2kT. Introducing this 
cutoff into the integral in Eq. (^) leads immediately to 
the logarithmic form for the fluctuation conductivity: 



2e2 7rK2T, kT 
a = — ; In rj/ K <C to. 
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(6) 



In the regime that m <^ "tj/k, the original integral in the 
conductivity is dominated by the minimum at x = 0. 
Performing the integrations in this regime results in the 
more standard form H for the conductivity: 



TO <C rj/n. 



(7) 



Should the transition to the superconducting state oc- 
cur at the finite temperature Tc and the effective field 
theory close to the transition is given by Eq. (0) with 
TO^ oa {T — Tc), we find that the fluctuation conductiv- 
ity is consistent with the standard Aslamazov-Larkin |^ 
form when the dissipation is of intermediate strength. 
However, in the weak dissipation limit, that is, weak in 
comparison to the inverse correlation length, the stan- 
dard 1/(T — Tc) is replaced by the logarithimic behavior 
sufficiently close to Tc. 

While the logarithmic conductivity follows naturally 
from an expansion in powers of r]/(Km) in the conductiv- 
ity, this appears to be the first time this result has been 
derived. Consider the regime to ^ kT. Fluctuations are 
suppressed in this regime. In this case, x± 3> 1 and its 
subsequent contribution to the integral is exponentially 
small. Thus, it will compete with the contribution from 
the vicinity of x = 0. Evaluating Eq. and retaining 
both contributions results in the total conductivity: 
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(8) 



The exponentially-decaying term is in agreement with the 
result of van Otterlo, et. al. ||2^, though their coefficient 
is a factor of 8 smaller. 

We now go beyond the quadratic theory we have con- 
structed. It is expedient to generalize to an N-component 
vector field, V'i'('^ik) and introduce the quartic term 
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V'l/ (wi , ki ) V^i. (w2 , k2 ) ('^a , ka ) V'a' (t^4 , k4 ) (9) 



into our free-energy functional. Because the phase order- 
ing transition is no longer determined by the vanishing 
of the coefficient of the quadratic term at fc = w = 0, it 
is expedient then to redefine 5 as the bare temperature 
independent coefficient of the quadratic term. 

We then introduce the auxiliary field, A(a;,k), to 
decouple the quartic term by means of a Hubbard- 
Stratonovich transformation. Formally, the auxiliary 
field, A(u;,k) is related to the correlation function, 
(tAqIw, k)'0* (w, k)). Consequently, if we invoke the 
saddle-point approximation, A(w,k) = \fP\5^.o5-k,Q, in 
the N ^ oo limit, the following self-consistency equation 
for the gap — S + \ where 



1 



(27r)2 m2 



(10) 



must be solved 



Let us define a = rj/'i-Kn'^T, h = 
rn? I^t^kT , and c = K/2'kk,T . Upon perform- 



ing the sum over boson frequencies and subsequent k- 
integration by introducing the relativistic hard cutoff, A 
in Eq. (foh, we obtain 



(11) 



n , C/T, T{a + ic)T{a-ic)K 
i"' = S In ■ 



2tt 



r{a + b)r{a-b)m 



as the self-consistency condition for the effective inverse 
correlation length. Here T{z) is the Gamma function. 
While Sachdev, Chubukov and Sokol have obtained a 
similar expression, the form derived here is more suitable 
for the asymptotic analysis of the three critical regimes 
of interest. 

Expansion of the gamma functions when T is the small- 
est parameter reveals that the behaviour of depends 
strongly on the sign of the quantity A = S + U{A + 
0{ri/ k)) / Attk. At the quantum critical point, A = and 
consequently m? = 0. At any finite temperature, the 
form of , which enters directly into the conductivity, is 
determined by the relative strength of three parameters: 
kT, a, and rj/ k. Consequently, the following regimes 
arise as depicted in Fig. (1). 
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FIG. 1. Heuristic phase diagram as a function of temper- 
ature T and A. The quantum critical point corresponds to 
A = 0. Regions I, II, and III correspond to 1) quantum dis- 
ordered, 2) quantum critical, and 3) renormalized classical. 
The dotted line demarcates the region in which dissipation 
dominates the physics of the correlation length. 



a) Quantum disordered (insulator): Two conditions 
(see Fig. (1)) demarcate this region: A > and 
A ^ kT. Expansion of Eq. ( pd] ) in this limit indi- 
cates for A rj/ K, m = AttkA/U and for rj/n ^ A, 
m = ^47r2?7A/([/ln {jp/nA)). We used the fact that 
from microscopic considerations of the model, k/U ~ 
0(1) close to the transition point . Substitution of these 
results into Eq. (^) leads to the corresponding conductiv- 
ity in the quantum disordered regime, b) Renormalized 
Classical: This regime is characterized by A < but 
|A| > kT as shown in Fig. (1). The solution to Eq. (|l|) 
in this limit results in the following forms for the gap 
parameter m: 




|A| > kT > t]/k 
|A|,77/k> kT 



(12) 



In the first sub-regime, m can be both smaller than or 
greater than 7]/n. Should m be greater than r^/n, then 
the logarithmic form for the conductivity applies. How- 
ever, because of the exponential dependence of to, the 
resultant conductivity 



4e2 



h r]U 



7]/ K < TO 



(13) 



is independent of temperature. Because rj/ k is small, the 
prefactor of e^ /h in the conductivity is large. This result 
is truly remarkable in light of the experiments by sev- 
eral groups ||l^,|l8|,^ that have reported a leveling-off of 
the conductivity in 2D thin film superconductors at low 
temperatures. Note that the conductivity as well as the 
width of the region over which <t is constant increase as 
the magnitude of the dissipation decreases. Hence, the 
model presented here is in principle compatible with the 
low values of the resistivity once the leveling-off ensues. 
Ultimately Mason and Kapitulnik |18| predict that dissi- 
pation leads to a "metallic" state at T = in field-tuned 
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IST's. However, within the framework presented here, a 
true "metahic" state does not exist at T = 0. We find, 
in fact, that at sufficiently low temperatures such that 



m < rj/ K, the second of Eqs. (12) apply and the conduc- 
tivity 



{eyh){T]/K^T) exp ( ^) kT > tj/h 
[e /n) exp I I ki <f]/K 



diverges exponentially as T ^ as is expected on the 
superconducting side. Nonetheless, the existing experi- 
ments cannot rule out the possibility that the resistivity 
ultimately vanishes at zero temperature. 

c) Quantum critical region. In this regime, kT ^ |A|. 
We consider for simplicity the transition point, A = 0. 
From Eq. (0), it follows that to = ^fl^^ fi{K^T /ri), 
where /i is a numerically calculable function. We find 
that in the limiting case nT <C rj/ m = \/2Tr ArjT, where 
A = 2.68723 is the larger of the roots of the equation 
r{x)\/x = ^/2tt. In the opposite limit, kT 3> rj/n, we 
recover the well-known jl^,^ result that to = QkT with 
Q = 21n((-\/5 -I- l)/2). In general, the conductivity is 
given by CT = (e^//i)/2(K^T/7y, k/J7), with /2 a universal 
function. To calculate /2, we investigate Eq. (^). For 
kT ^ ti/k, the main contribution to the integral over x 
comes from the vicinity of x± leading to a conductivity 
of the form, a = 2Ki{e^ /h){K^T/f]) with 



t^dt 



(i2 + 62/4) sinh^(i2 + 02/4) 



0.0149. (15) 



For kT <C ?7/k, it is convenient to perform the k- 
integration in Eq. (^) first. We obtain for the conduc- 
tivity, a = {2K2/n){e'^ /h) where 



K. = 



dt 



sinh t 



IT A ( TT -kA 

1 arctan 

t \2 t 



= 0.0076. 



We see that in the limit T 0, the conductivity acquires 
a universal value independent of the dissipation. The in- 
dependence of the critical resistance on 77 is a consequence 
of hyperuniversality |^0|. One should note, though, that 
the extremely small numerical prefactor, is likely to be 
an artifact of this model. 

Our analysis of the role of dissipation in the vicinity 
of the quantum critical point has revealed a key new fea- 
ture, namely the leveling of the resistivity at sufficiently 
low temperatures. To determine if this leveling is exper- 
imentally relevant, we evaluate the resistivity by invert- 
ing Eq. (|l^. As shown previously U/k = 7/2 and 
K — l/4:Ec = 0.3K^^ close to the quantum critical 
point. Here Ec is the charging energy. Consequently, 
the level resistivity reduces to i? = Ro{7 /2TT){ri/ k)/\\A\, 
where Ro = h/Ae^ — QAk^l. Typical values for the lev- 
eled resistance are R = lOVl |jl^ which would require 
the ratio r?/K|A| = .01. For |A| in the range [0.01,.l], 



we find that the corresponding value of the dissipation 
ranges between 10~^ < rj/ n < 10~^ which is certainly in 
the weak dissipation regime as required by the present 
analysis. From the condition to « rj/ n, we find from Eq. 
( p!^ ) that the leveling behavior ceases at T O{10mK). 
Hence, we find that for realistic parameters, our analysis 
predicts a level resistance that is consistent with experi- 
mental trends. 
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